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ABSTRACT 

We report a dynamically evolving low ionization broad absorption line flow in the QSO SDSS 
J133356. 02+001229.1 (at z e m ~ 0.9197). These observations are part of our ongoing moni- 
toring of low ionization broad absorption line (BAL) QSOs with the 2m telescope at IUCAA 
Girawali observatory (IGO). The broad Mg n absorption with an ejection velocity of 1.7xl0 4 
km s _1 , found in the Sloan Digital Sky Survey (SDSS) spectra, has disappeared completely 
in our IGO spectra. We found an emerging new component at an ejection velocity of 2.8xl0 4 
km s" 1 . During our monitoring period this component has shown strong evolution both in its 
velocity width and optical depth and nearly disappeared in our latest observations. Acceler- 
ation of a low velocity component seen in SDSS spectrum to a higher velocity is unlikely 
as the Mg n column densities are always observed to be higher for the new component. We 
argue that the observed variations may not be related to ionization changes and are consistent 
with absorption produced by multi-streaming flow transiting across our line of sight. We find 
a possible connection between flux variation of the QSO and iV(Mg n) of the newly emerged 
component. This could mean the ejection being triggered by changes in the accretion disk or 
dust reddening due to the outflowing gas. 
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1 INTRODUCTION 

It is widely believed that the kinetic energy output from quasars 
through the outflows are as important as their radiative output. They 
are believed to play an i mportant role in regulating the g rowth of 
supermassive blackholes (Silk &Reesll 19981: iKh ig 2003) and star 
formation in the host galaxies (e.g., iBower et alJl2006T) . In addi- 
tion, outflows influence the enrichment of the surrounding inter- 
galactic medium. Therefore, it is important to understand the ori- 
gin and evolution of QSO outflows. Large scale outflows manifest 
itself in the form of broad absorption lines (BALs) in QSO spec- 
tra. BAL QSOs comprise of up to 40% of the total QSO population 
jDai et all 2008, and references therein). The correct explanation 
for the observed incidence of BAL QSOs is s till a subject of debate 
between the orientational (see for example, [Hines & Wills! fl995l : 
iGoodrich & Millerlll995tlMurrav et all 1995l;lElvisl2qO0h and evo- 
lution ary models (see forexample. lHazard et aljl 984; Beck er et al.l 
l200d) . 

Line variability studies of BAL QSOs are useful for under- 
standing the physical conditions and dynamics of the gas close to 
the central engine. The time variability of C iv and Si iv absorp- 
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tionis reported in several cases (for example, ISrianand & Petitieanl 
1200 it iLundgren et"al]|2007l ; iGibson et alj|2008h . Such a variabil- 
ity could be related either to variations in the ionization state or to 
the covering factor of the absorbing gas. However, the most inter- 
esting cases are the ones where the flow emerges afresh or shows 
strong dynamical evolution (i.e variation in the absorption profile 
and signatures of acceleration). There are five previous reports of 
emerging C iv BAL discovery in the liter ature [TEX 1726+34 4 
Jm3 l2002h , SPSS J105400.40+ 034801.2 faamaim et ail 120081) , 
WPVS 007 ( lLeighlvetalJl2009h. Ton 34 jKrongold et al.l l201Ch 
and PG0935+417 I Rodrigu ezHidalgo et alj|20 111) 1. These studies 
attribute the observed dynamical evolution to multiple streaming 
wind moving across the line of sight. Detecting an emerging Mg n 
flow will be very interesting as Mg n BALs are considered to be a 
possible man ifestation of a Q SO's efforts to expel a thick shroud of 
gas and dust JVoit et all 1993b . 

Here, we report the first discovery of an emerging 
Mg ii outflow in a QSO SDSS J133356.02 +001229.1(hereafter 
11333+0012) at a redshift z cm ~ 0.9197 faewett et al] 1200 it 
iTrump et al.l2006h . This source is part of our sample of low ioniza- 
tion BAL QSOs that are being spectroscopically monitored at IU- 
CAA Girawali observatory (IGO). Our spectra revealed the emer- 
gence of a new Mg n broad absorption component that was not 
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Table 1. Log of Observations 



Source 


Instrument 


Date 


MJD 


Exposure Time 


A Coverage 


Resolution 


S/N" 










(min) 


(A) 


(kms" 1 ) 






SDSS 


28-04-2000 


51662 


45x1 


3800-9200 


150 


27 




SDSS 


15-02-2001 


51955 


105x1 


3800-9200 


150 


40 




IGO/IFOSC 7 


03-04-2008 


54559 


40x2 


3800-6840 


300 


11 




IGO/IFORS 1 


26-02-2009 


54888 


45x3 


3270-6160 


360 


21 




IGO/IFOSC 7 


26-03-2009 


54916 


45x3 


3800-6840 


300 


12 


J1333+0012 


IGO/IFOSC 7 


19-04-2009 


54940 


45x2 


3800-6840 


300 


5 




IGO/IFORS 1 


22-01-2010 


55218 


45x1 


3270-6160 


360 


19 




IGO/IFORS 1 


25-01-2010 


55221 


45x1 


3270-6160 


360 


16 




IGO/IFORS 1 


14-03-2010 


55269 


45x3 


3270-6160 


360 


16 




IGO/IFOSC 7 


06-04-2011 


55657 


45x3 


3800-6840 


300 


19 



" calculated over the wavelength range 5800-6200A 



present in the Sloan Digital Sky Survey(SDSS) spectra obtained 7 
years earlier. We have been monitoring this new Mg n absorption 
component for the past 4 years. In section 2, we provide details 
of the observations and data reduction. Results and discussions are 
presented in Sections 3 and 4 respectively. 



2 OBSERVATION AND DATA REDUCTION 

Our spectroscopic observations were carried out with IUCAA Faint 
Object Spectrograph (IFOSC) mounted on a 2m telescope at IU- 
CAA Girawali observatory (IGO). The details of these observa- 
tions together with that of the archival SDSS spectra are given 
in Table [U We have used Grism #1 and Grism #7 of IFOSC in 
combination with a long slit having a width of 1.5 arcsec. This 
combinations have the wavelength coverage of 3270-6160 A and 
3800-6840A respectively. Typically the observations were splitted 
in to exposures of 45 minutes. Spectra obtained on closely sep- 
arated dates were combined for better signal to noise. Cleaning 
of the raw frames and ID spectral extraction were carried out us- 
ing the standards procedures in IRArfJ We opted for the variance- 
weighted extraction with "doslit" procedure. Wavelength calibra- 
tions were done using standard helium neon lamp spectra and flux 
calibrations were done using a standard star spectra observed on the 
same night. Air-to-vacuum conversion was applied before adding 
the spectra. Individual spectra were combined using 1/cr 2 weight- 
ing in each pixel after scaling the overall individual spectra within 
a sliding window. The error spectrum was computed taking into 
account proper error propagation during the combining process. 

Fig.Q]shows the SDSS spectrum taken on 15-02-2001. The 
expected positions of broad emission lines (at z em =0.9197) are 
marked with blue dotted vertical lines. Narrow Mg n, Mg i and 
Fe ii absorption lines associated with the intervening absorbers at 
z a bs~ 0.898 and z a b s ~ 0.836 are distinctly detected. In addition, two 
broad Mg n absorption components (marked in Fig.[T]as 'B' and 'R' 
to denote the blue and red components respectively) are clearly de- 
tected. We do not detect any associated Fe n lines and the Al in lines 
are redshifted out of the cov erage. A visual compari son of this spec- 
trum with the Fe template of Vestergaard & Wilkes (2001) suggests 
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Figure 1. The SDSS spectrum of J1333+0012. The expected positions of 
broad emission lines are marked with blue dotted lines. The narrow ab- 
sorption lines at z a bs~ 0.898 and 0.836 are shown by black and red dashed 
lines. Arrows show the red (marked as R) and blue (marked as B) com- 
ponents of Mg n BAL discus sed in this article. Fe emission template from 
IVestergaard & Wifkesl I2OO1I) is shown in the upper panel for comparison. 

that the continuum of J1333+0012 may have considerable contri- 
bution from broad Fe emission lines. 



3 RESULTS 

iTrump et all d2006|) have fitted the spectrum of J1333+0012, ob- 
served on 15-02-2001, using a template spectrum. Their best fitted 
template resulted in a spectral index ( f v oc V') of a v = 1.64 and 
E(B-V) = 0.295 for a SMC type dust. They have identified only one 
distinct BAL component at z a bs= 0.8 1 1 (corresponding to R compo- 
nent defined above). They assigned a Balnicity index and Absorp- 
tion index of 36±0.42 and 459±0.75 respectively to this source. 
Apart from this R component (that has consi stent profile between 
two SDSS spectra and the KECK spectra of iHewett etall J200lh 
observed on 30-04-2000) we find a new component at z a b s = 0.764 
(defined above as B) in the 2001 SDSS data (not clearly detected 
in the previous epoch SDSS spectrum) having a velocity width of 
-6000 km s" 1 . 

In Fig. [2] the IGO spectra obtained at different epochs are 
compared with a reference SDSS spectrum obtained in year 2000 
(i.e MJD 51662). For displaying purpose we have scaled the fluxes 
to match the flux of the reference spectrum. In the bottom half of 
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Figure 2. The spectra of J1333+0012 obtained at different epochs are compared with the initial SDSS spectrum. The observed flux is scaled to match the SDSS 
flux. The vertical lines show the wavelength range of R and B components as in Fig.[T] The ratio between the SDSS and the corresponding IGO spectrum is 
plotted in the lower panel in each plot. A lower order polynomial fit to this curve and the associated error are also shown. 



each panel we show the ratio of spectrum obtained at a given epoch 
to that of the reference spectrum. This curve is smooth apart from 
the wavelength range of R and B components. We normalise this 
curve with a smooth lower order polynomial (shown as red dotted 
curve in Fig. [2) to remove any large wavelength scale continuum 
flux differences between the two spectra. 

In our IGO 2008 (MJD 54559) spectrum (see panel b in Fig©, 
the R component ( that was seen with near identical optical depth in 
two SDSS spectra) has completely disappeared. The R component 



has never reappeared again in any of our subsequent IGO spectra 
(see Fig[2j- At the same time, the B component is seen to have 
varied in depth and widened in velocity in the 2008 (MJD 54559) 
data. In 2009 (MJD 54888), the B component further deepened in 
the absorption troughs. But, no changes in the velocity structure are 
seen (see panels c and d in Fig.©. It diminished in strength in 2010 
and has almost vanished in 201 1 (see panels e and f in Fig. [2}. 

We calculated the Mg n column density fo r both R and B com- 
ponents using apparent optical depth method dSavage & Sembachl 
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Table 2. Column Density Measurements 



MJD 


R Component 




B Component 






log N(Mg ii ) J 


Av° 


z b 

abs 


log N(Mg ii) 1 


Av" 


z b 

abs 




(cm' 2 ) 


(kms- 1 ) 




(cm" 2 ) 


(kms" 1 ) 




51662 


13.62+ 0.20 


2829 


0.811 


<13.42* 


7000 




51955 


13.75+0.13 


3934 


0.812 


13.36±0.25 


6143 


0.764 


54559 


<13.57* 


3000 




14.26±0.39 


12303 


0.747 


54888 


<13.20* 


3000 




14.50±0.10 


12909 


0.746 


54916 


<13.42* 


3000 




14.46±0.20 


11900 


0.745 


55218 


<13.09* 


3000 




14.28±0.13 


11474 


0.754 


55657 


<13.35* 


3000 




13.76±0.36 


10016 


0.760 



The error in A^Mg n ) does not include the contribution from continuum 
placement uncertainties; " v elocity range that co ntains 90% of the 
integrated optical depth fsee lLedoux et alj2006h : b optical depth weighted 
redshift; * 3 <x Upper Limit of A^Mg n ). 




52000 



53000 54000 
MJD 



55000 



Figure 3. Top panel: Light curve (in Johnson's V magnitude) of 
J1333+0012 from CRTS and SDSS (first two points) measurements. Black 
points are the average of closely spaced observations and the associated er- 
ror is obtained from errors in individual measurements. The x-error bars 
corresponds to the time period over which the individual measurements are 
averaged. The epochs of spectroscopic observations are marked as blue dot- 
ted vertical lines. Bottom panel: the measured W(Mg n ) of the B component 
vs MJD. 



Il99lh assuming the gas to be optically thin and covering the back- 
ground source completely. In case of the component B, we have 
used the normalised ratio spectra discussed in Fig. [2] to estimate 
iV(Mg ii ). For the component R that is relatively narrow we use the 
observed flux in the neighborhood of the absorption line for contin- 
uum normalization. The measured Mg n column densities (or 3tr 
upper limits), velocity width (Av) and z a bs for both the components 
are summarized in Table [2] 

In Fig.[3]we plot the photometric light curves of J1333+0012 
obtained in the Catalina Real-Time Transient Survey (CRTS, 
iDrake et alj2009l) . These observations were carried out with an un- 
filtered set up and were converted to V magnitude using 10 to 100 
well selected G-type dwarf c alibratio n stars found in each frame 
and the method described in Henden (2000). The zero point un- 
certainty is found to be < 0.08 mag. The two left most points in 
the light curve corresponds to the transformed SDSS photometric 



and fiber magnitu des respectivel y. We use the transformation equa- 
tion given in jjester et alll2005h for z < 2. 1 QSOs for this. Fiber 
magnitude should be considered as an upper limit on the actual 
magnitude. The QSO has dimmed between the early epoch CRTS 
observations and IGO 2008 observations, reached a minimum in 
2009 and then brightened thereafter. The maximum seen amplitude 
of the variation is ~0.2 mag. Very interestingly, a similar but oppo- 
site variation is seen in the iV(Mg n ) of the B component during 
this period, i.e., when the quasar flux was minimum, the measured 
N(Mg ii) is maximum. 



4 DISCUSSION 

The main result of our monitoring of J1333+0012 is that the R 
component of Mg n absorption has disappeared completely after 
the year 2001 and the B component emerged, became stronger and 
wider and subsequently disappeared between the years 2001 and 
2011. There is also tentative evidence for a possible connection 
between the QSO flux and /V(Mg n) measured for the B component. 
Now we discuss various scenarios to explain these observations. 

It is very tempting to say that the R component got acceler- 
ated to the B component and later vanished from our line of sight. 
From Table [2] it is clear that the value of N(Mg n) and Av of com- 
ponent B in later epochs are much higher than what is seen in the 
R component in earlier epochs. If we use thin shell approximation, 
then mass conservation suggests that ^A^H) should be conserved, 
where r and N(H) are the distance of the shell from the QSO and 
total hydrogen column density respectively. We expect r to increase 
with time so A'(H) (and hence Af(Mg n) should decrease with time). 
The observed A'(Mg n) values do not not support this scenario. 

One could argue that there were always two fluid components 
and what we see is due to the variations in the ionization parameter. 
Above we have seen that the QSO has varied by ~ 0.2 mag. How- 
ever the required chang e in A^Mg n) is 0.6 to 1 dex. It is clear from 
Fig. 2 of lHamannl d 19971) that required change in A'(Mg n) can not 
be produced by 0.2 dex change in ionization parameter. In addition 
the change in the Mg n optical depth of R and B components are in 
the opposite directions. This scenario can be consistent if the vari- 
ation in UV continuum is much larger than the 0.2 mag seen in the 
V band and the initial ionization parameters of the two components 
are very different. We do not favor this scenario based on the lack 
of large colour variations between our IFOSC spectra and the near 
consistency of Mg n emission line. 

In all the previous case s of emerging flows 



Jm3 l2002l: iHamannetalJ 120081: iKrongold et alj |2010|; 
iRodriguez Hidalgo et alj|201 luHall et alJuOl lh the most favoured 
explanation is the transverse motion of the absorbing gas. In 
particular the evolution of the B component can be explained as the 
absorbing gas transiting through our line of sight as the component 
has emerged, evolved and disappeared. To consider this possibility 
we derive some basic parameters for J1333+0012. 

We obtained the B ba nd magnitud e of this Q SO using the mea- 
sured u and g mag nitudes Ijester et alj (s ee l2005i) . We use this and 
the prescription of iMarconi et alj 1 20041) . to get the bolometric lu- 
minosity, L fco/ = 9.39 x 10 46 ergs s -1 . The corresponding mass ac- 
cretion rate, M h , is, 



M, 



Lboi 
~7a 



l6.9M a /yr 



(1) 



where e, the mass to energy conversion efficiency is taken as 0.1. 
The blackhole mass derived fro m Eddington accretion is 8xl0 8 
M Q . Following lHall et al.1 JioTTh . we find the diameter of the disc 
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within which 90% of the 2700-A continuum is emitted D2700 ~ 
10 16 cm. 

The variability of the R component has occurred over a period 
3994 days in the observers frame or 208 1 days in the rest frame of 
the QSO. We can use this as an upper limit on the transit time of the 
absorbing gas across the continuum emitting region. If we assume 
the projected transverse size of the absorbing gas is much smaller 
than the emitting region (assumed to be a face-on disk), then we 
can estimate the transverse velocity v± > 550 km s -1 . We find the 
radial velocity up to v r = 2.8xl0 4 km s~'and the radial displace- 
ment of the gas during the transit is ~ O.lpc. Even if the projected 
transverse size of the gas is 10 times the continuum emitting re- 
gion, the expected transverse velocity will be much lower than the 
measured radial velocity. 

The main question is, is the emergence of new component in 
any way related to the structural changes in the accretion disk?. 
The possible connection between QSO dimming and N(Mg n) is 
interesting. It could either mean that the ejection being triggered 
by some events in the accretion disk that caused reduction in the 
accr etion efficiency or extinction by dust in the transiting gas. 

iTrump et al. (2006) identified this system as BAL based on the 
presence of Mg n absorption in the R component. The latest spec- 
trum of J1333+0012 taken with IGO is nearly devoid of any broad 
Mg n absorption. We do not detect Mg i, Fe n and Ca n absorp- 
tion associated with either of the broad Mg n components. It will 
be interesting to study the properties of other strong high ionization 
lines using UV spectroscopy. Continuous optical monitoring of this 
source to see recurrence of new outflow episode will be useful for 
modeling this interesting source. 



5 CONCLUSIONS 

In this paper, we report a dynamically evolving Mg n broad ab- 
sorption line outflow in SDSS J133356.02+001229.1. The R com- 
ponent detected in the SDSS spectra taken in 2000 and 2001 was 
never seen in the IGO spectra taken in the latter epochs. The blue 
component B emerged in 2001. It continued to increase its optical 
depth till 2009 and started to diminish thereafter. This component 
has almost vanished in the most recent spectra in 201 1. 

We provide arguments for not favoring the acceleration of a 
single cloud to a higher velocity and changes in the ionization con- 
ditions causing the observed spectral variability. The observed vari- 
ability can be best explained by multiple streaming gas moving 
across our line of sight. An outflowing gas component with trans- 
verse velocity one tenth of the observed radial velocity and flow 
extent up to 10 times the size of the accretion disk can easily ex- 
plain these observations. 

We find the QSO has dimmed by about 0.2 mag (in V-band) 
during the monitoring period. This may be the manifestation of the 
emerged outflow being triggered by the structural changes in the 
accretion disk or due to reddening by dust in the absorbing gas. 
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